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SYNTHETICALLY FOCUSED 
ULTRASONIC DIAGNOSTIC IMAGING SYSTEM 
FOR TISSUE AND FLOW IMAGING 

5 This application is a continuation in part of US 

patent application serial number 10/136,880 of the 
same title, which was filed on April 30, 2002. 

This invention relates to ultrasonic diagnostic 
imaging systems and, in particular, to synthetically 
10 focused ultrasonic diagnostic imaging systems which 

image tissue, flow, and flow spectra. 

Today's conventional ultrasound systems produce 
images by scanning a region of the body with 
ultrasound beams electronically formed by a 
15 transducer array. The timing of signals applied to 

the array elements determines the direction in which 
a beam is transmitted and the depth of focus of a 
beam. A plurality of adjacent beams are transmitted 
to adequately spatially sample the target region. 
20 The echoes from along the beam directions are 

processed in known ways to form an image of the plane 
or volume scanned by the beams and depict 
characteristics of the tissue or bloodflow within the 
plane or volume . 
25 when a beam is transmitted, the beam profile 

determined by the transducer element timing can be 
set to produce the best focus around a selected focal 
_.- - j~„-iv^ ^or-Nt-Vi a transmitted beam can 
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have only one such focal point; once the beam is 
30 launched, the focal characteristic cannot be changed 

or extended. The clinician will generally set the 
focus to the depth at which the anatomy of interest 
may be found. If the clinician wants the optimal 
focus to be at two or more depths or range of depths, 
35 multizone focus must be used, whereby two or more 
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beams are transmitted in the same direction, each 
focused at a different depth. The entire image plane 
or volume is scanned by multiple beams in each 
direction, and the best focused portions of the sets 
5 of beams are spliced together to produce a composite 

image which exhibits the optimal focus at a plurality 
of depths. However, each differently focused beam 
set will extend the time required to acquire the 
image information, which reduces the frame rate of 
10 display. Going from one focal zone to two can 

approximately halve the frame rate, for example. 

During reception, optimal focus can be obtained 
dynamically over the full imaging depth. This is 
because the receive focus can be dynamically adjusted 
15 electronically by continually adjusting delays in the 

receive beamformer as echoes are received over the 
full range of depth. The transmitted beams are 
constrained to a single focal point or region since 
the transmit focus is accomplished acoustically and 
20 not electronically. 

A technique for overcoming this transmit focus 
limitation is known as synthetic focus. Synthetic 
focus is described in US Pat. 4,604,697, for 
instance. In a synthetic focus system, each 
25 transducer element or subset of transducer elements 

is actuated sequentially. The transmission from each 

element or group of elements is not focused 

_ __..„*-■! ~,n v K,it- pnvprs i-hfi entire imaqe region. The 
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echoes from each transmission are received by all of 
30 the elements concurrently and stored. These echoes 

are then combined in different combinations with 
different effective delays, thereby forming coherent 
echoes at points in the image region which are 
effectively focused at all points. This technique 
35 produces image data which is in effect optimally 
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focused on both transmit and receive, providing 
images which are optimally focused throughout the 
entire depth of field. In such an approach it is not 
necessary to first beamform onto distinct A- lines, 
5 although this can be done. Rather, the echo data can 

be directly reconstructed onto image pixels. 

While the l 697 patent illustrates a synthetic 
focus technique which is limited to imaging 
stationary tissues, it is also desirable to extend 
10 the technique to flow imaging. US Pat. 5,349,960 

illustrates one extension of a synthetic focus 
architecture to Doppler. In this patent a synthetic 
focus acquisition subsystem is used to produce real 
and imaginary parts of a complex signal, which is 
15 then Doppler processed. The full synthetic focus 

data set is not used. Rather, reception only occurs 
on the transmitting element. It is desirable to have 
a synthetic focus architecture which utilizes more 
than one receive element and which produces data for 
20 all ultrasonic modes, including B mode, colorflow, 

and spectral flow, without the need for hardware or 
software dedicated to only a specific mode of 
operation . 

The parent application addresses these desires 
25 by providing an ultrasonic diagnostic imaging system 

in which a one or two dimensional transducer array is 
operated to acquire synthetic focus ultrasound 

. . , -, _ t . ,n„„ ^r,« rvr a subset of the elements 
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of the array and receiving echo signals on all of the 
30 array elements in response to a transmit event. The 

received signals are digitized and stored in a 
memory. Sets of delays are applied to combinations 
of the stored signals to produce motion maps for a 
plurality of different velocities. One or more of 
35 the motion maps are used to produce B mode, 
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colorflow, spectral, and other ultrasound image modes 
which cannot be produced by conventional ultrasound 
systems. The images, either two or three 
dimensional, are in optimal focus throughout the 
5 image plane or volume and can be produced at high 

frame rates of display. Various transducer element 
sequencing may be employed to increase the frame 
rates even higher. 

In accordance with the principles of the present 

10 invention, the synthetic focused ultrasound system is 

operated in a hybrid fashion in which operation 
alternates between synthetic focus acquisition and 
conventional focused beam acquisition. This makes 
possible, for example, the acquisition and display of 

15 harmonic images for contrast and/ or tissue harmonic 

imaging. Speckle artifacts in the synthetic focused 
ultrasound images may be reduced by combining signals 
from different look directions. For instance, sets 
of motion maps may be produced for different sub- 

2 0 apertures of the array transducer, then compounded to 

- reduce speckle. One or more identified regions of 
interest within a synthetic focused ultrasound image 
may be processed differently from other regions of 
the image to highlight or better define particular 
25 motional characteristics within the regions of 

interest, such as turbulent flow or different 
velocities of flow or motion. 
In the drawings : 

FIGURE 1 illustrates in block diagram form the 

3 0 major functional components of a typical ultrasound 

system of the prior art; 

FIGURE 2 illustrates the presence of two 
stationary targets in the field of a synthetically 
focused transducer array; 
35 FIGURE 3 illustrates the presence of one 
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stationary and one moving target in the field of a 
synthetically focused transducer array; 

FIGURES 4a-4e illustrate the raw data acquired 
by scanning the two stationary targets of FIGURE 2; 
5 FIGURES 5a-5f illustrate synthetic focus 

processing of the raw data of FIGURES 4a-4e with 
delays chosen for the condition of no motion at the 
location of one of the targets; 

FIGURES 6a-6f illustrate synthetic focus 
.0 processing of the raw data of FIGURES 4a-4e with 

delays chosen for the condition of motion at the 
location of one of the targets; 

FIGURES 7a-7e illustrate the raw data acquired 
by scanning the two targets of FIGURE 3; 
L5 FIGURES 8a-8f illustrate synthetic focus 

processing of the raw data of FIGURES 7a- 7e with 
delays chosen for the condition of no motion at the 
location of one of the targets; 

FIGURES 9a-9f illustrate synthetic focus 
20 processing of the raw data of FIGURES 7a-7e with 

delays chosen for the condition of motion at the 
location of one of the targets; 

FIGURE 10 illustrates in block diagram form a 
synthetic focus data acquisition system suitable for 
25 acquiring the data illustrated in FIGURE 4-9; 

FIGURE 11 illustrates an ultrasound system 
constructed in accordance with the principles of the 

^v-oonnl- -J rnran t" i OT1 • 

FIGURE 12 illustrates the operation of a linear 
3 0 array transducer to acquire synthetic focus 

ultrasound data; 

FIGURE 13 illustrates a second technique for 
operating a linear array transducer to acquire 
synthetic focus ultrasound data; 
35 FIGURE 14 illustrates the operation of a two 
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dimensional array transducer to acquire synthetic 
focus ultrasound data for a high frame rate of 
display with reduced temporal artifacts; 

FIGURE 15 illustrates the time interleaving of 
5 synthetic focus data acquisition and conventional 

focused beam data acquisition in a hybrid embodiment 
of the present invention; 

FIGURE 16 illustrates the use of sub-apertures 
to acquire synthetic focused data sets for speckle 
10 reduction by spatial compounding; and 

FIGURE 17 illustrates a synthetic focus 
ultrasound image in which different processing is 
performed for synthetic focus data that is applicable 
to a region of interest in the image. 
15 Current diagnostic ultrasound systems may be 

described as having three distinct operating modes : 
B-mode echo (i.e., grey-scale), spectral Doppler 
(i.e., single-gate Doppler with high spectral 
resolution in only a single sample volume) , and color 
20 Doppler (i.e., multi-gate Doppler where a single 

scalar quantity such as mean Doppler frequency, or 
total Doppler power, is estimated in many sample 
volumes but with essentially no spectral resolution.) 
This partitioning is reflected in both the system 
25 architecture and system functionality. Typically, 

specialized hardware is required to implement each 
mode and only one mode is operational at once. An 
example of this standard ultrasound system 
architecture is shown in FIGURE 1. A probe 10 has an 
3 0 array transducer 11 which transmits beams along 

selected lines through the image field under control 
of a beamformer 12. Echoes from along the lines are 
formed into lines of coherent echo signals by the 
beamformer 12. The echoes are applied to a signal 
35 processor 14 which has dedicated hardware and 

ATL-34 0 



software for the various modes of operation. A B 
mode processor 14a will process the echoes for grey- 
scale display or a Doppler processor 14b will process 
the echoes for colorflow or spectral Doppler display. 
Often separate processors are used for colorflow and 
spectral Doppler processing. The processed echo 
information then undergoes image processing 16, and 
the resultant image is displayed on an image display 
18. So-called "simultaneous" operation actually 
involves rapid temporal multiplexing between 
different ones of these modes. Such partitioning is 
necessary in conventional systems in order to 
optimize performance in each mode while maintaining 
acceptable frame rates. For example a typical 
spectral Doppler processor would operate too slowly 
to process all of the pixels for a real time 
colorflow display. 

However, such partitioning is neither 
fundamentally necessary nor desirable. Ideally, an 
ultrasound system would simultaneously estimate the 
amount of backscattered power in all velocity 
increments or "bins" at all points in the image, and 
at frame rates sufficient to track temporal changes. 
From this general solution for the properties of the 
target field, conventional operating modes could be 
derived. For example, standard B-mode images are 
obtained from the amount of power in the zero 

v>-; ^ pa! nnnnl imaaes mav be computed 
by integrating over non-zero velocity bins to 
estimate mean velocity or total power, spectral 
Doppler displays would show the discrete values in 
all bins, and so on. It is this approach which is 
taken with the present invention. 

FIGURES 2-9 illustrate examples of synthetic 
focusing, which is utilized to acquire raw echo 



340 



-8- 



information for processing in accordance with the 
principles of the present invention. FIGURE 2 
illustrates a transducer array 11 which, for ease of 
illustration, is shown with only five transducer 
5 elements 1-5. The center axis 6 of the image field 

extends from the central element 3. Located in the 
center of the image field at spatial location (43,0) 
is a stationary target 7. Located to the left of the 
center axis 6 as viewed from the array 11 is another 
10 stationary target 8 at spatial location (45,10). 

A synthetic focus data set is acquired as shown 
in FIGURE 4. Each set of lines of echoes 4a, 4b, 4c, 
4d, and 4e illustrates the echoes acquired by all of 
the transducer elements 1-5 in response to 
15 transmission by a single element. The set of lines 

4a shows echoes received from targets 7 and 8 in 
FIGURE 2 in response to transmission by element 1. 
The leftmost echo on each line is the echo from the 
central target7 . The symmetrical distribution of 
20 this group of echoes 20 reveals that the target 7 is 

in the center of the image field, and its earlier 
reception in time reveals that it is at a nearer 
spatial location (with a shorter round-trip echo 
path) than the target 8. The asymmetrical time 
25 distribution of the second (later in time) group of 

echoes reveals that target 8 is off -axis and at a 
greater depth than target 7. The earlier occurrence 
of echo 22 at element 1 than that of echo 24 at 
element 5 reveals that target 8 is closer to element 
30 1 than element 5. That is, target 8 is on the left 

side of the image field. This fact is also revealed 
by the later occurrence of the second group of echoes 
as successive elements 2-5 are used as the transmit 
element in echo sets 4b, 4c, 4d, and 4d, as the 
35 round-trip distance of each echo path increases as 
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the transmitting element is further removed from 
target 8 with each successive transmit element 2-5. 

The round- trip distance from any transmit 
element to any receive element for any location in 
the image field is a priori known. Given that the 
speed of sound in the medium is known everywhere or, 
as is commonly the case, can be approximated by a 
constant, sets of delays can be applied to the raw 
echo data shown in FIGURE 4 to reveal the foregoing 
characteristics of the targets in the image field. 
Moreover, the delays can be further applied for 
different assumptions of motion by the targets. 
FIGURE 5 shows the results when a set of delays are 
applied for location (45,10) in the image field for 
the assumption of a static target, that is, it is 
assumed that target 8 is not moving during the 
transmit sequence. When these delays are applied to 
the echo data, it is seen that all of the echoes from 
the different transmits line up and can be coherently 
summed, as shown by coherent echo 26' on line 5f. 
These delays, which are appropriate to align echoes 
from a static target at location (45,10), are 
inappropriate for a target at (43,0), and hence the 
noncoherence of the echoes in group 20 from target 7 
is increased. The summation of these echoes will 
produce only a minor ripple on line 5f as can be seen 
on line 5f below echo group 20. The amplitudes of 

■, j - ~„ i-ino arp normalized bv the 
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number of signals used in the summation (25 in this 
example.) The waveform of 5f is then detected in a 
known manner for B-mode echoes to reveal the 
scattering strength of target 8 . 

FIGURE 6 illustrates the application of a second 
delay set to the raw echo data, a delay set which is 
appropriate for the assumption that, during the 
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transmit sequence, target 8 moves from location 
(45,10) to location (49,0) as shown in FIGURE 3. 
When each successive element in the array is used as 
the transmit element, this delay set assumes, the 
target 8 is at a different location on its way to 
location (49,0) as indicated by the open circles in 
FIGURE 3. However, target 8 is not moving in this 
case. The result of this delay set is that coherence 
will be produced for the set of echo lines 6a, for at 
this time the target 8 is assumed to be at its 
starting location (45,10). However, the delay sets 
for the echo line set 6b, 6c, 6d, and 6e each assume 
that the target 8 has moved to a location which is 
incrementally closer to final location (49,0). In 
fact, the target 8 is static and has not moved. As a 
result, the second group of echoes in each line set 
becomes progressively more uncorrelated as FIGURE 6 
illustrates. It is also seen that the first set of 
echoes 20 in echo line set 6e are almost in time 
coincidence, since target 7 is on the central axis 
and the delay set used for line set 6e assumes that 
target 8 has moved to an on-axis location but at a 
slightly greater depth than that of target 7. Since 
this assumption, that target 8 is moving, is 
incorrect, the summation line 6f reveals no coherent 
echoes for the delay set for this assumption. 

FIGURE 7 shows a raw echo data set for the 
condition where target 7 is stationary and target 8 
is moving from location (45,10) to location (49,0) 
during the transmit sequence interval shown in FIGURE 
3. As was the case for FIGURE 4, echo line set 7a 
depicts echoes received by all elements 1-5 (a-e) 
when element 1 is the transmit element; echo line set 
7b depicts echoes received by all elements when 
element 2 is the transmit element; echo line set 7c 
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depicts echoes received by all elements when element 
3 is the transmit element; and so forth. At the time 
of each transmit event the target 8 is at one of the 
successive circled locations including or between 
locations (45,10) and (49,0), as shown in FIGURE 3. 

FIGURE 8 illustrates the results obtained when a 
delay set is applied to the raw data which is 
appropriate for a static target at location (45,10), 
the same delay set that was applied in FIGURE 5. The 
second echoes from target 8 will be in time coherence 
as shown in echo line set 8a for this delay set, 
since at this time the target 8 is at the assumed 
location (45,10). But for successive echo line sets 
8b- 8e the second group of echoes 26 become 
increasingly noncoherent, for at each successive 
transmit time the target 8 is at a greater distance 
from the assumed location. The echoes from the 
static target 7 always demonstrate noncoherence , as 
was the case in FIGURE 5. Consequently the summation 
line 8f shows no coherent echo signals for this delay 
set and assumption of no motion. 

Finally, FIGURE 9 illustrates the result 
obtained when the delay set chosen is that which 
matches or tracks the motion of target 8 as it moves 
from location (45,10) to location (43,0) over the 
transmit interval of five pulses. For that set of 
delays all of the second group of echoes 26 from 

, 4. o -,n -i^r, ^r>v.a->-<=nt-l v as shown in FIGURE 9, and 
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by the coherent echo on the summed line 9f. The 
echoes 20 from target 7 remain noncoherent and do not 
produce significant signal content on line 9f, except 
from echo line set 9e where both targets are on the 

center line 6. 

As can be seen from the above, a delay set is 
effectively the coefficient set for a matched filter 



340 



-12- 



which will reveal a characteristic of the image field 
when objects in the image field have characteristics 
matching that of the coefficient set. When a target 
moves in both the direction and at the speed for 
which the delay set is effective, that motion is 
revealed. Furthermore, it is seen that a delay set 
has a unique motional characteristic in both speed 
(velocity) and direction. Thus, matched filters, 
each having a characteristic velocity and direction, 
a vectorial characteristic, can be used to detect 
motion in the image field with those vectorial 
characteristics. In the example shown in FIGURE 9, 
the direction of motion is neither purely axial nor 
purely lateral. Rather, the velocity has an oblique 
characteristic . 

An ultrasound system acquisition subsystem for 
acquiring echo line sets as illustrated above is 
shown in FIGURE 10. A transducer array 11 of M 
elements has M channels 1-M connected to the elements 
of the array. In the illustrated embodiment each 
channel includes a transmitter as shown by TX1 in 
channel 1 for the first element. However, since only 
one element or subgroup of elements are activated at 
any one time, one or a few transmitters could be 
multiplexed to each element when that element is to 
transmit. The transmit pulse from the transmitter is 
coupled to the element which is to transmit by way of 

■ / „^ ^,.t-;*-^v. 70 wh-irh orotects the 
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receive circuitry from energy from the transmitter. 
The echo signals received by the element following 
transmission are coupled by the transmit/receive 
switch to an amplifier 34 which may be gain 
controlled to apply time gain compensation to the 
received echo signals. The echo signals are 
digitized by an A/D converter and may then optionally 
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be demodulated into in-phase and quadrature signal 
components by a quadrature bandpass filter (QBP) 38, 
which may also impose a desired passband on the 
received signals. The QBP is described in more 
detail in US Pat. 6,050,942. 

The echo line sets following transmission by an 
element of the array are stored in a memory 40. For 
the case of a one dimensional array used in 2D 
(planar) imaging, the memory has three dimensions: 
m, n, and time. In the illustrated embodiment each 
column has M rows, one for each element that receives 
echoes following a transmit event. The memory has N 
columns, one for each of N transmit events. 
Dimension n can also be thought of as a "slow time" 
dimension. If each element is used to transmit, N 
will equal M. However, it is possible to transmit 
with fewer than all elements, which will provide a 
reduced signal to noise ratio and, depending on the 
choice of elements, reduced lateral spatial 
resolution but a higher frame rate, as the full data 
set used for image production can be acquired in less 
time. It is not necessary to actuate the transmit 
elements in any particular sequence. Elements at 
separated distances along the array can be actuated 
in succession, skipping back and forth along the 
array. Once the predetermined set of echo line sets 
has been acquired, the data is used for image 

1 .i.J — « .J ~ ^. ,~i -i V>qH hoi r\\kT 
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Consider the case of a conventional Doppler 
system which uses a 64 element phased array 
transducer (M = N = 64) with a center frequency fo = 
3 MHz. The velocity resolution is related to the 
frequency resolution Af according to the Doppler 
equation: 

Av = (c Af) / (2 fo) 
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where Af is determined by the reciprocal of the FFT 
window length (assuming a rectangular window), i.e.: 
1/(L AT), where L is the number of samples (transmit 
pulses) in the FFT window, and AT is the period 
between pulses. A typical value might be AT = 
200uSec which produces a pulse repetition frequency 
(PRF) of 5 kHz. A typical value for L is 64 so that 
the conventional Doppler system has a velocity 

resolution of 

c / (2 fo L AT) =2.0 cm/sec. 
Now consider operation of the same 64 element 
transducer in a synthetic focus system in which each 
element is pulsed once. The total observation 
interval is N AT. Further assume that the motion is 
axial. Then, if the difference, Av, between the 
velocity that a matched filter is tracking and the 
actual target velocity is such that the target moves 
an additional X/2 during the observation interval, 
then the echoes from the target undergo a full cycle 
of phase shift and the output of the matched filter 
will be minimized. Therefore, Av represents the 
velocity resolution and equals: 
Av = a/2)/ (N AT) 
= c/ (2 fo N AT) 
25 =2.0 cm/Sec. 

This is the same as for conventional Doppler. Thus 
there need be no loss in velocity resolution in a 
synthetic focus approach. 

The above discussion considers axial velocity 
resolution. Lateral velocity resolution will 
generally not be as great since the rate of change of 
phase with position is typically less in the lateral 
plane. Lateral velocity resolution can be increased 
by using larger transmit and receive apertures. The 
dynamic transmit focus possible in synthetic focus 
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allows the use of large transmit apertures without 
depth of focus problems. 

Of course, the velocity resolution of the 
synthetic focus system can be increased merely by 
extending the observation interval. This can be 
achieved by increasing AT or, preferably, by pulsing 
each element multiple times during the observation 
interval. Although this reduces frame rates, the 
result is only to reduce a very high frame rate to a 
fully acceptable frame rate. A conventional B-mode 
image typically uses on the order of 100 transmits 
along 100 distinct lines. For the value of AT used 
in this example, each frame would require 20 mSec (or 
a frame rate of 50 Hz) . In colorflow mode, an 
15 additional 3 0 Doppler lines might be used. Each 

Doppler line would typically require perhaps 10 
transmits increasing the total acquisition time by 60 
msec and reducing the overall frame rate to 12.5 Hz. 
in contrast, for this example, the frame rate for 
synthetic focus would be set by N = 64 transmits and 
would be 78.1 Hz regardless of whether the output 
mode was B-mode or colorflow. 

On the other hand, if reduced velocity 
resolution is acceptable then fewer pulses (and 
therefore fewer transmit elements) can be used. This 
results in even higher frame rates, which makes the 
synthetic focus approach one worth considering for 
real time 3D imaging. 

In an analogous fashion to spectral analysis, 
apodization functions can be applied across the 
observation windows in order to attenuate sidelobes 
in the velocity response. If the acquisition 
sequence is monotonic with respect to the transmit 
element position (that is, the sequencing steps 
across from one side of the transmit aperture to the 
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other) , then the apodization function serves double 
duty - it also attenuates sidelobes of the spatial 
(beam) response. 

Various sources of "spectral broadening" may 
5 reduce velocity resolution beyond that predicted by 

the narrow band analysis above. A major cause of 
spectral broadening in conventional Doppler is 
"transit time broadening". With short pulses or 
narrow beams (large apertures) , the targets only 
10 remain in the sample volume for a limited period, 

thus reducing the observation interval. One solution 
is to lengthen the pulse and/or broaden the beam 
(reduce the aperture) but this reduces spatial 
resolution. The present invention avoids this 
15 problem since it inherently "tracks" the targets, 

that is, the sample volume moves with the targets. 
Therefore, it is possible to achieve both high 
spatial resolution and greatly reduced spectral 
broadening . 

20 Remaining sources of spectral broadening include 

transit time in the elevation dimension (when using 
one -dimensional arrays) and, to a lesser extent, the 
secondary effects of changes in the RF point spread 
function with position. However, since the targets 

25 will likely be tracked over relatively small 

distances, this latter effect is not a major problem. 
While "wall filtering" is critical in color 

_ -, ,j 4-^,/= r^^i- wei nri t-v resolution of the 
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short ensemble and because color Doppler measurements 
30 are mean velocity estimates) , its function in 

spectral Doppler (which has much better velocity 
resolution) is largely to reduce down stream dynamic 
range requirements and allow the use of FFT-based 
spectral analysis. It is therefore the case that, 
35 provided care is taken to control velocity sidelobe 
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response and to have adequate processing dynamic 
range, wall filters are not fundamentally necessary 
to implement the present invention. However, if some 
degree of wall filtering is desired, it may be 
implemented simply by transmitting from each element 
more than once. Then, the wall signal can be 
estimated and subtracted from the return from each 
individual transmit element before application of the 
delay sets. The velocity response of this process 
can be tailored through suitable choice of the 
intervals between transmit events and filter weights 
and can be different for different motion maps. For 
example, no wall filtering would be used when 
calculating the maps corresponding to static or 
slowly moving targets. 

In an embodiment of the present invention 
spatial resolution can be at least as good as that of 
conventional B-mode and Doppler imaging because: 1) 
Synthetic focus allows dynamic transmit focus; 2) It 
is not necessary to use narrow bandwidth bursts (poor 
axial resolution) to reduce spectral broadening; and 
3) It is not necessary to use wide beams (poor 
lateral resolution) to reduce spectral broadening. 

Frame rates are determined by the number of 
transmit pulses. Synthetic focus allows the entire 
frame (not just one line) to be reconstructed from 
these pulses. Therefore, very high frame rates can 

, . j mv,^ -i ~ r^-f nart-inilar sianificance for 
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color flow imaging and, as mentioned above, for 3D 
imaging from 2D arrays. It is emphasized that it is 
not necessary to compromise other aspects of system 
performance (such as field of view or spatial 
resolution) to achieve high frame rates. Only 
velocity resolution has a direct effect on frame 
rates but, as the example above shows, it is possible 
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to achieve much higher frame rates than conventional 
color flow while at the same time achieving much 
greater velocity resolution. 

Sensitivity is a major issue in Doppler studies 
since scattering from erythrocytes is weak. 
Sensitivity is fundamentally limited by the amount of 
energy that is transmitted and the amount of thermal 
noise contributed by the transducer elements and 
receive amplifiers. This places ' synthetic focus at a 
disadvantage by approximately the square root of M 
(in terms of voltage signal to noise ratio) if the 
voltages applied to the transmit elements are not 
changed relative to conventional transmit 
beamforming. However, when it is recalled that 
synthetic focus transmits are non-focused and have 
low per-element duty cycles it is apparent that the 
voltage applied to each transmit element can be 
significantly increased beyond that used for the 
element in conventional beamforming transmission 
before whatever limiting factor (e.g.,: SPTA, max MI, 
scanhead temperature, saturation, etc.) is reached. 
Analysis based on peak pressure levels then predicts 
a square root of M advantage in sensitivity for 
synthetic focus over conventional beamforming. Coded 
waveforms, with large time-bandwidth products, may 
also be employed to increase sensitivity. 

in the embodiment shown in FIGURE 10, the raw 
echo data from synthetic focus acquisitions as 
illustrated in FIGURE 10 is processed by K delay sets 
(coefficient sets) to form K motion maps, each for a 
different motion characteristic. Selected ones of 
these motion maps are then used individually or in 
combination to produce all of the standard ultrasound 
imaging modes, and, in addition, if desired, modes 
not possible in conventional systems. 
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Mathematically, this process begins with the 
acquisition of a set of waveforms f n , m (t) , where 
indices n and m identify a particular 
transmit/receive element pair in the full data set. 
A set of K images of the form g k [u) is constructed, 
where u is a position vector in the ultrasound image 
u = (x,y) for a 2D image and u = (x,y,z) for a 3D 
image) and k specifies the k th motional 
characteristic. The K images are processed using a 
motion analyzer which extracts the type of motion 
and/or image format which is to be displayed. An 
individual image is then displayed, or a number of 
images are shown simultaneously for a compound 
display format such as colorflow. An entirely new 
set of waveforms may be acquired to produce the next 
image in the real time sequence, or a partial set 
acquired and processed with previously acquired 
waveforms to produce images at an increased frame 
rate but with some degree of temporal persistence 
from frame to frame. In the latter embodiment, the 
motion maps 60 of FIGURE 11 must be stored in 
undetected format. However, memory 4 0 can be reduced 
in size or even eliminated. 

The images in this process may be expressed as: 

N M 

where W n>in is an apodization weighting function and 

. , . , . _i_ j — ~ -p-,,-^ -i of motion 
T nmk is a delay weignc wmcn x B a j.wc~w- — 

and', for a constant velocity may be expressed as 

pj r „ S-(n-\)PRIv k \+fr*m -u-(n-\)PRJv k \ 



T..> =• 



L n,m,k q 

A Tn is the location of the transmit element for the 
n th ' transmit event and A R , m is the location of the m fc 
receive element. PRI is the pulse repetition _ 
interval and the velocity of the motion map is v 4 



,th 
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If the velocity is not constant but accelerative, 
then the linear position term (n-l)PRIv, must be 
replaced by a more general function a (n,k) to 
describe the motion. 

One problem with the simple expression for the 
T nmk delay weight presented above is that the targets 
move slightly between the time that the n th transmit 
event occurs and the time that the target actually 
intercepts and scatters the energy from that 
transmit. This results in a slight focusing error 
(because the net displacement of the target from its 
assumed position is different for different transmit 
and receive element pairs) . The error becomes more 
significant with increasing target velocities. A 
more precise calculation of the delays would take 
this effect into account. 

A functional block diagram of an ultrasound 
system constructed in accordance with the principles 
of the present invention is shown in FIGURE 11. The 
elements of a transducer array, which may be a ID 
array or a 2D array, are coupled to M 
transmit/receive channels (see FIGURE 10) which cause 
individual ones of the elements to transmit for 
reception by a selected group of receive elements. 
Typically, all of the array elements will receive 
echoes in response to each transmit event. The 
transmit/receive channels are coupled to the memory 
40 where the sets of echo waveforms are stored as 
described above. The stored echo information is 
operated upon by K delay or coefficient sets 50 and 
detected in the conventional manner to produce K 
motion maps 60, each for a different velocity vector. 
A motion analyzer 7 0 then selects various ones of the 
motion maps to produce the desired mode or modes of 
ultrasound display. For instance, the motion map for 
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v t =0 can be used to produce a grey-scale (B-mode) 
image 84. However, when motion is present, contrast 
will be reduced identically to conventional synthetic 
focus imaging. It is preferred that the B-mode image 
84 be constructed from either the peak pixel values 
over all motion maps or, preferably, the integration 
of the magnitude or power in each pixel over all 
motion maps in order to avoid reduction in contrast 
due to motion. A B-mode tissue image is then 
overlaid with a colorized motion image made by 
appropriately integrating the non-zero velocity 
images above a certain velocity determined by the 
desired wall filter cut-off to compute mean velocity 
(normalized first moment), the result being a 
colorf low image 86. The velocity data at any pixel 
or group of pixels can be plotted graphically to 
produce a spectral flow display 82, or a group of 
such displays 96 at selected points in the image. An 
embodiment of the present invention can produce the 
equivalent of a multigate spectral Doppler display at 
any or all points in the image from a single 
acquisition sequence. The velocity values can be 
accumulated from the non-zero motion maps to produce 
a total power display 95, the equivalent of a power 
Doppler image. By spatially displaying velocity 
information of tissue using a wall filter setting 
which passes low velocities (with the possible 

„ f--iacno mntion imaqe 92 
exception ot zero veiuciLj,, ~ — 

can be produced. A power motion image 91 (see US Pat. 
5,718,229 entitled "MEDICAL ULTRASONIC POWER MOTION 
IMAGING") can be produced from static motion maps 
from successive acquisitions. Since the motion 
vector at every pixel is present in the motion 
images, a vector flow image 93 can be produced. An 
image 94 showing the flow velocity with the greatest 



340 



-22- 



power can be produced if desired. Other display 
formats will readily occur to those skilled m the 
art . 

in a constructed embodiment the system 
resolution will be determined by the number of 
velocity bins (for discrete motion maps) at each 
pixel in a given direction and the number of 
directions for which velocity bins are produced. For 
example in the axial direction there may be ±64 bins 
for velocities both toward and away from the 
transducer, for a total of 128 bins. In the lateral 
direction where resolution is inherently poorer there 
ma y be only ±16 or 32 bins. The number of directions 
for which velocities are estimated between these 
orthogonal directions is also a variable to be 
chosen, and can be four, eight, or sixteen 
directions, depending upon the desired directional 
resolution. In a two-dimensional imaging system in 
which the vector velocity is solved at every point m 
the image, 32,000 or more bins per pixel may be 
needed. A three dimensional system will increase 
this to millions of bins. 

As mentioned above, the acquisition rate can be 
increased by transmitting with less than all of the 
elements. For instance, FIGURE 12 illustrates a ID 
array HI in which only the end elements 112 and 114 
are used for transmission. For each transmit event, 

f^-rablv used 
all of the elements or tne duay *. . 

for reception, although fewer may be employed. 
Additional elements in the center of the array may be 
used for transmission as desired, or any other sparse 

array combination. 

The acquisition sequence embodied by FIGURE 12 
can have the same high spatial resolution as a 
sequence which uses all elements to transmit . 
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However, when different targets have different 
velocities there may be undesirable cross-talk, or 
ambiguity, between the spatial response and the 
velocity response. In other words, the ability of 
the system to distinguish between a first target or 
group of targets at a first position moving at a 
first velocity and a second target or group of 
targets at a second position moving at a second 
velocity is diminished. 

This ambiguity can be controlled by the choice 
of which transmit and receive elements are used and 
in which order. For example, the ambiguity can be 
minimized by using an acquisition sequence of the 
type illustrated in FIGURE 13 . Here an inner pair of 
elements 102 and 104 are used to sequentially 
transmit. Element group 106 is used to receive when 
the transmit is from 102 and element group 108 is 
used to receive when the transmit is from 104. This 
sequence has the property that the mid-points between 
transmit element 102 and all the elements in receive 
aperture 106 are co-located with all the mid-points 
between transmit element 104 and all the elements in 
receive aperture 108. Other sequences that preserve 
the commonality of the transmit/receive mid-points 
are also possible. For example, the separation 
between the pair of transmit elements can be 
increased (or reduced) while the number of elements 

i -; vn^npoH ( n-r increased) 

in the receive aptJiuuic j.s — — 

appropriately. 

Sequences of the type illustrated in FIGURE 13 
reduce the ambiguity between range and velocity at 
the expense of reducing spatial resolution whereas 
the sequence in FIGURE 12 maximizes spatial 
resolution at the expense of increased ambiguity. 
Other tradeoffs, where there is an intermediate 
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degree of mid-point commonality, are also possible. 
For example, the transmit element pair separation of 
FIGURE 13 could be used with the full receive 
apertures of FIGURE 12. Or, the transmit sequence 
can use more than two elements. 

in general, provided that the received data is 
preserved in memory, the present invention allows 
such a tradeoff to be made after data acquisition has 
occurred. The tradeoff can be very flexible. For 
example, different motion maps can be constructed 
using different numbers and locations of transmit and 
receive elements. It is generally preferred to use 
more elements, and hence higher spatial resolutions 
for low velocities. For higher velocities such as 
those encountered when imaging relatively weak 
scatterers such as blood, it is generally desirable 
to reduce cross-talk from slower moving but stronger 
scattering targets such as tissue. Here it is 
preferred to trade off some spatial resolution for 
reduced cross -talk by using a sequence which 
increases commonality of the mid-points. Such 
tradeoffs can be made adaptively, based on analysis 
of the spatial and velocity characteristics of the 
particular target situation. In addition, the 
results from multiple, different, acquisition 
sequences for a given velocity can be detected then 
combined non-linearly . Such combinations can also be 
adaptive . 

FIGURE 14 illustrates a 2D array for three 
dimensional imaging. Only the elements around the 
periphery of the array are shown in detail; the 
elements internal to the periphery of the 2D array 
are represented by the crossed dots in the center of 
the drawing. Transmission from this 2D array 121 for 
synthetic focus imaging may proceed by transmitting 
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from ones of the peripheral elements in succession 
and receiving by all of the array elements following 
each transmit event. For instance, corner element 
122 may be activated first, followed by transmission 
from the other corner elements 124, 126, and 128 in 
succession. These four transmissions provide 
sufficient data to create motion maps with full 
spatial resolution. However, all receive elements, 
which typically number in the thousands, must be 
used. Alternatively, if the progression of transmit 
elements then proceeds around the periphery by 

-i ^ m =,»-,)-a 1^? 134 136, and 138 next, and 
actuating elements li^, ±Ji0 ' 

so on until each of the peripheral elements has 
produced a transmit event, then reception need only 
occur on the elements on the periphery of the array 
to create full spatial resolution motion maps. The 
sequence is then repeated for the next full echo data 
set The frame rate may be the inverse of the time 
required to acquire a full data set from transmission 
by all of the peripheral elements, or may be 
increased by updating the image when part of a new 
sequence has been acquired. 

Since the synthetic focus echo information is 
buffered, that is, it is stored in memory, the same 
data that is used to derive motion information and 
images may also be used to estimate and correct for 
tissue aberrations. Techniques for estimating and 

former delay 
correcting tor aDeriat-iuiio ~- - 

control are described in US patents 6,023,977, 

6,131,458, 6,027,447 and 6,223,599 for example. 

However, in all of these embodiments estimates for 

aberration correction are produced from a current 

data set, then used during the reception of 

subsequent data sets. In a synthetic focus system of 

the present invention, since the transmit beamformmg 
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is done electronically rather than acoustically, 
transmit corrections may be applied to the same data 
set used to estimate the corrections. This avoids 
the need for real-time estimation and re- transmission 
5 and also can provide more accurate aberration 

correction. 

FIGURE 15 illustrates an embodiment of the 
present invention in which hybrid acquisition is 
employed. That is, the ultrasound system alternates 
10 between acquiring synthetic focus echo data and 

standard focused beam echo acquisition. Such a 
system makes possible the production of harmonic 
images (tissue or contrast) in a synthetic focus 
setting. There are two limitations to harmonic 
imaging in a conventional synthetic focus system. 
First, the intensity of harmonic signals, 
particularly tissue harmonic signals, is very low 
relative to fundamental frequency signals. This is 
because harmonic generation is proportional to the 
square of the acoustic pressure at the target. Since 
synthetic focus generally only uses a few elements of 
a transducer array for each transmission, and may 
only use a single element, the pressure of the 
acoustic energy field thus developed is far below 
that which is present at the focus of a standard 
focused transmit beam, and may be insufficient to 
produce diagnostically useful harmonic images. 

. •,. - u;>^« r\-f ar*Yir> sianals etmoloved 

Furtnermore, me Cuuiuj.nin a ~. 

in the production of a synthetic focus image is 

generally done linearly and is predicated on assumed 

linearity in the acoustic field. While this 

assumption is true for fundamental signals it is not 

the case for nonlinear higher harmonics. 

Hybrid acquisition as shown in FIGURE 15 

addresses these problems and enables harmonic imaging 
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in the context of synthetic focus imaging. In this 
embodiment synthetic focus signal acquisition is 
performed during a first time interval from t„ to t x 
These signals may be acquired and buffered until all 
the data needed for synthetic focus motion imaging 
has been acquired. The synthetic focus data from the 
transducer elements is then combined to form a motion 
or flow map such as described above. This can be the 
flow image component of a color flow image, for 
example. During a subsequent time interval from , t 4 to 
t 2 conventional focused beams are transmitted and the 
echoes acquired at a harmonic frequency. This 
harmonic echo data may be used to produce the tissue 
harmonic image component of a color flow image The 
flow and tissue components are overlaid or displayed 
together as a colorflow image in which the tissue is 
imaged harmonically. In the third interval from t 2 to 
t, the system performs synthetic focus signal 
acquisition again, followed by focused beam 
transmission and acquisition during the fourth 
interval from t, to t 4 . The third and fourth 
intervals may comprise the acquisition time for a 
second color flow image frame, or a number of time 
interleaved intervals greater than two may be 
alternated to acquire both synthetic focus and 
focused beam signals for a display image. A hybrid 
embodiment such as this may, under certain 

conditions, be faster tnan 4 w n««n~ - - 

system; can do motion tracking as described above; 
and enables images to be reprocessed from buffered 
data for different sets of parameters such as 
beamforming characteristics, wall filtering, or 

thresholding. 

A synthetic focus motion imaging system may also 
be used to produce quantified data. Quantification 
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may be applied to the synthetic transmit aperture 
images. In particular, parameters that are not 
routinely available, such as motion direction 
(vector) , linear and angular accelerations can be 
computed from synthetic focus data. 

in a synthetic focus system it is not necessary 
to process all of the data from the entire aperture 
coherently. For example, coherent processing may be 
performed on multiple sub-apertures and the results 
combined incoherently to reduce speckle through 
spatial compounding. With buffered synthetic focus 
data sets, this processing can be applied and the 
degree of coherence adjusted for processing or 
reprocessing of the same acquired data set. In the 
embodiment of FIGURE 16, the full aperture of a 
transducer array 11 is divided into three sub- 
apertures identified as sub-aperture 1, sub-aperture 
2, and sub-aperture 3. The signals received by each 
sub-aperture are synthetically focused as described 
above to produce a set of one or more motion image 
maps for each sub-aperture. The corresponding motion 
image maps from the different sub-apertures are then 
compounded together to reduce speckle by virtue of 
the different look directions of each sub-aperture. 
If desired, an image map produced from the full 
aperture 11 can be compared or blended with the 
compounded image map set of the sub-apertures to 

. __t j — ^^finaiiv useful imaae data, 
produce aaaiuiuiidi mayn^^v^^j 

Further aspects of spatial compounding may be found 
in US patents 6,210,328, 6,126,599, and 6,135,956. 

It is not necessary to process the entire 
frequency spectrum coherently, either. If coherent 
(pre-detection) processing is performed on multiple 
sub-bands of the received synthetic focus echo 
signals and the results combined incoherently, 



speckle can be reduced by frequency compounding. As 
in the previous embodiment, with buffered synthetic 
focus data sets this processing can be applied and 
the degree of coherence adjusted for processing or 
reprocessing of the same acquired data set. Further 
aspects of frequency compounding are described in US 

patent 4,561,019. 

In a synthetic focus motion imaging system, 
different regions of the image, termed regions of 
interest or ROIs, can be processed differently 
according to analysis of local characteristics. The 
variability can be controlled by a user interface 
which allows the user to position and size a ROI box. 
FIGURE 17 illustrates an embodiment of the present 
invention in which the synthetic focus data pertinent 
to a ROI is processed differently from data pertinent 
to the surrounding image field. FIGURE 17 shows a 
synthetic focus ultrasound image 140 displaying an 
image of the heart 142. The heart has an upper 
chamber 146 and a lower chamber 144 connected by a 
heart valve 148. Because the blood flow in the 
vicinity of the heart valve 148 may be uniquely 
different from that of the heart chambers, an ROI box 
150 is positioned around the heart valve 148. For 
instance the blood flow around the valve may be more 
turbulent than the more laminar flow of the 
surrounding chambers and blood vessels. The blood 

. ■ i... _ji 4-v,« v=iTro ids mav exhibit 
flow in tne vicmxoy ui >-"^ - 

different velocities than the surrounding blood flow. 
Accordingly, the synthetically focused data pertinent 
to the ROI box 150 is differently processed in 
consideration of these characteristics. The data 
from a different number of transducer elements may be 
used to synthetically focus the motional 
characteristics within the ROI box, or different 
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spatial resolution may be employed. If the 
velocities within the ROI box are anticipated to be 
higher than the surrounding velocities, lower 
resolution may be employed within the ROI box to 
5 follow or track the heart valve blood flow more 

accurately. The buffered synthetic focus data 
pertinent to the ROI box may be reprocessed to reveal 
particular characteristics such as different linear 
velocities and/or acceleration. The set of motion 
10 maps used to produce the image within the ROI box may 

be different than the set of motion maps used outside 
of the ROI box. Thus, synthetic focusing is 
performed differently within and outside of the ROI 
box in consideration of the diagnostic conditions 

15 encountered. 

Synthetic focus motion data can also be 
subjected to post detection image processing such as 
XRES, an image processing technique available from 
Philips Ultrasound, which provides a more smoothly 
enhanced tissue texture in an ultrasound image. 

Synthetic focus motion imaging may also benefit 
from wider aspect ratio and flat panel displays such 
as described in US patent 6,516,215. Since synthetic 
focus motion imaging can produce simultaneous images 
25 in different display formats or motion 

characteristics from the same synthetic focus data 
set, displays such as these, with greater display 

— oimHipHr focus imaqes to be 

areas, penuiL uiuxux^j-v- ^j^^ -- 

displayed at the same time for diagnoses based on 
multiple diagnostic parameters or the comparison of 
multiple images. 

In summary, an embodiment of the present 
invention can provide one or more of the following 
advantages : 
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• Good spatial resolution of echoes in all 
modes, with improvement due to dynamic transmit 
focusing; 

High frame rates, even in color flow mode; 
Reduced spectral broadening; 
High velocity resolution; 

Improved (or reduced requirement for) wall 
iltering for colorflow imaging; 
New display modes for flow; 

Estimation of vector velocity with no angle 
correction, isotropic color displays; 

Harmonic as well as fundamental imaging; 
Speckle reduction; 

Images with spatially variable processing; 



or 



A unified system architecture with no 



partitioning for each mode 
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